Studying the Brain Monoaminergic Systems and Neurotrophic Factors in Minipigs with High and Low Tolerance to the Presence of Human.

Here, we present the first evidence for brain adaptation in pigs tolerant to the human presence, as a behavioral trait favoring domestication. The study was carried out on minipiglets from population bred at the Institute of Cytology and Genetics (Novosibirsk, Russia). We compared the behavior, metabolism of monoaminergic neurotransmitter systems, and functional activity of the hypothalamic-pituitary-adrenal system, as well as neurotrophic markers in the brain of minipigs differing by tolerance to human presence (HT and LT - high and low tolerance). The piglets did not differ in the levels of activity in the open field test. However, the concentration of cortisol plasma was significantly higher in minipigs with a low tolerance to the presence of humans. Moreover, LT minipigs demonstrated a decreased level of serotonin in the hypothalamus and augmented levels of serotonin and its metabolite 5-HIAA in the substantia nigra as compared to HT animals. In addition, LT minipigs showed increased content of dopamine and its metabolite DOPAC in the substantia nigra and decreased dopamine level in the striatum as well as reduced content of noradrenaline in the hippocampus. Increased mRNA levels of two markers of the serotonin system - TPH2 and HTR7 genes - in the raphe nuclei and in the prefrontal cortex, respectively, were associated in minipigs with a low tolerance to human presence. However, the expression of genes regulating a dopaminergic system (COMT, DRD1, and DRD2) in HT and LT animal groups varied depending on brain structure. In addition, a decrease in the expression of genes encoding BDNF (brain-derived neurotrophic factor) and GDNF (glial cell line-derived neurotrophic factor) was revealed in LT minipigs. The results may contribute to our understanding of the initial stage of domestication in pigs.

Pathogen exposure influences immune parameters around weaning in pigs reared in commercial farms.

BACKGROUND: Multiple antigenic stimulations are crucial to immune system training during early post-natal life. These stimulations can be either due to commensals, which accounts for the acquisition and maintenance of tolerance, or to pathogens, which triggers immunity. In pig, only few works previously explored the influence of natural exposition to pathogens upon immune competence. We propose herein the results of a multicentric, field study, conducted on 265 piglets exposed to contrasted pathogen levels in their living environment. Piglets were housed in 15 different commercial farms, sorted in two groups, low (HSLOW)- and high (HSHIGH)-health status farms, depending on their recurrent exposition to five common swine pathogens. RESULTS: Using animal-based measures, we compared the immune competence and growth performances of HSLOW and HSHIGH pigs around weaning. As expected, we observed a rise in the number of circulating leucocytes with age, which affected different cell populations. Monocyte, antigen-experienced and cytotoxic lymphocyte subpopulation counts were higher in piglets reared in HSLOW farms as compared to their HSHIGH homologs. Also, the age-dependent evolution in γδ T cell and neutrophil counts was significantly affected by the health status. With age, circulating IFNα level decreased and IgM level increased while being greater in HSLOW piglets at any time. After weaning, LPS-stimulated blood cells derived from HSLOW piglets were more prone to secrete IL-8 than those derived from HSHIGH pigs did. Monocytes and granulocytes issued from HSLOW pigs also exhibited comparable phagocytosis capacity. Altogether our data emphasize the more robust immunophenotype of HSLOW piglets. Finally, piglets raised under higher pathogen pressure grew less than HSHIGH piglets did and exhibited a different metabolic profile. The higher cost of the immune responses associated with the low farm health status may account for lower HSLOW piglet performances. CONCLUSIONS: Altogether, our data, obtained in field conditions, provide evidence that early exposure to pathogens shapes the immune competence of piglets. They also document the negative impact of an overstimulation of the immune system on piglets' growth.

On an association between fear-induced aggression and striatal-enriched protein tyrosine phosphatase (STEP) in the brain of Norway rats.

Striatal-enriched protein tyrosine phosphatase (STEP) is a signal transduction protein involved in the pathogenesis of neuropathologies. A STEP inhibitor (TC-2153) has antipsychotic and antidepressant effects. Here, we evaluated the role of STEP in fear-induced aggression using Norway rats selectively bred for 90 generations for either high aggression toward humans (aggressive rats) or its absence (tame rats). We studied the effects of acute administration of TC-2153 on behavior and STEP expression in the brain of these animals and the influence of chronic treatment with TC-2153 on the behavior and STEP expression in aggressive rats in comparison with classic antidepressant fluoxetine, which is known to exert antiaggressive action. Acute TC-2153 administration decreased the aggressive reaction to humans in aggressive rats, while having no impact on the friendly behavior of tame rats. Moreover, in the elevated plus-maze test, the drug had an anxiolytic effect on both aggressive and tame rats. Aggressive rats demonstrated elevated levels of a STEP isoform (STEP46) as compared to tame animals, whereas acute TC-2153 administration significantly reduced STEP46 protein concentration in the brain of aggressive rats. Chronic treatment of aggressive rats with either TC-2153 or fluoxetine attenuated fear-induced aggression. Chronic administration of fluoxetine enhanced the exploratory activity in the elevated plus-maze test and decreased the STEP46 protein level in aggressive rats' hippocampus, whereas chronic TC-2153 administration did not affect these parameters. Thus, STEP46 can play an important role in the mechanisms of aggression and may mediate antiaggressive effects of TC-2153 and fluoxetine.

Genetic Factors Mediate the Impact of Chronic Stress and Subsequent Response to Novel Acute Stress.

Individual differences in physiological and biobehavioral adaptation to chronic stress are important predictors of health and fitness; genetic differences play an important role in this adaptation. To identify these differences we measured the biometric, neuroendocrine, and transcriptional response to stress among inbred mouse strains with varying degrees of genetic similarity, C57BL/6J (B), C57BL/6NJ (N), and DBA/2J (D). The B and D strains are highly genetically diverse whereas the B and N substrains are highly similar. Strain differences in hypothalamic-pituitary-adrenal (HPA) axis cross-sensitization were determined by plasma corticosterone (CORT) levels and hippocampal gene expression following 7-weeks of chronic mild stress (CMS) or normal housing (NH) and subsequent exposure to novel acute restraint. Fecal CORT metabolites and body and organ weights were also measured. All strains exposed to CMS had reduced heart weights, whereas body weight gain was attenuated only in B and N strains. Acute stress alone produced larger plasma CORT responses in the D and N strains compared to the B strain. CMS paired with acute stress produced cross-sensitization of the CORT response in the N strain. The N strain also had the largest number of hippocampal transcripts with up-regulated expression in response to stress. In contrast, the D strain had the largest number of transcripts with down-regulated expression following CMS and acute stress. In summary, we observed differential responses to CMS at both the physiological and molecular level among genetically diverse strains, indicating that genetic factors drive individual differences in experience-dependent regulation of the stress response.

Impact of Genetic Variation on Stress-Related Ethanol Consumption.

BACKGROUND: The effect of stress on alcohol consumption in humans is highly variable, and the underlying processes are not yet understood. Attempts to model a positive relationship between stress and increased ethanol (EtOH) consumption in animals have been only modestly successful. Our hypothesis is that individual differences in stress effects on EtOH consumption are mediated by genetics. METHODS: We measured alcohol consumption, using the drinking-in-the-dark (DID) paradigm in females from 2 inbred mouse strains, C57BL/6J (B6) and DBA/2J (D2), and 35 of their inbred progeny (the BXD family). A control group was maintained in normal housing and a stress group was exposed to chronic mild stress (CMS), consisting of unpredictable stressors over 7 weeks. These included predator, social, and environmental perturbations. Alcohol intake was measured over 16 weeks in both groups during baseline (preceding 5-week period), CMS (intervening 7-week period), and post-CMS (final 4-week period). RESULTS: We detected a strong effect of CMS on alcohol intake. A few strains demonstrated CMS-related increased alcohol consumption; however, most showed decreased intake. We identified 1 nearly significant quantitative trait locus on chromosome 5 that contains the neuronal nitric oxide synthase gene (Nos1). The expression of Nos1 is frequently changed following alcohol exposure, and variants in this gene segregating among the BXD population may modulate alcohol intake in response to stress. CONCLUSIONS: The results we present here represent the first study to combine chronic stress and alcohol consumption in a genetic reference population of mice. Differences in susceptibility to the effects of stressful environments vis-à-vis alcohol use disorders would suggest that the differences have at least some basis in genetic constitution. We have also nominated a likely candidate gene underlying the large individual differences in effects of stress on alcohol consumption.

Effects of divergent selection upon adrenocortical activity on immune traits in pig.

BACKGROUND: The sustainability of farming and animal welfare requires the reconsideration of current selection schemes. In particular, implementation of new selection criteria related to animal health and welfare should help to produce more robust animals and to reduce anti-microbial use. The hypothalamo-pituitary-adrenocortical (HPA) axis plays a major role in metabolic regulation and adaptation processes and its activity is strongly influenced by genetic factors. A positive association between HPA axis activity and robustness was recently described. To explore whether selecting pigs upon HPA axis activity could increase their robustness, a divergent selection experiment was carried out in the Large White pig breed. This allowed the generation of low (HPAlo) and high (HPAhi) responders to adrenocorticotropic hormone administration. RESULTS: In this study, we compared 23 hematologic and immune parameters of 6-week-old, HPAlo and HPAhi piglets and analysed their response to a low dose of lipopolysaccharide (LPS) two weeks later. At six weeks of age, HPAhi piglets displayed greater red blood cell and leucocyte number including CD8α+ γδ cells, cytotoxic T lymphocytes, naive T helper (Th) cells and B lymphocytes as compared to HPAlo individuals. The ability of blood cells to secrete TNFα in response to LPS ex vivo was higher for HPAhi pigs. At week eight, the inflammatory response to the LPS in vivo challenge was poorly affected by the HPA axis activity. CONCLUSIONS: Divergent selection upon HPA axis activity modulated hematologic and immune parameters in 6-week-old pigs, which may confer an advantage to HPAhi pigs at weaning. However, HPAlo and HPAhi piglets did not exhibit major differences in the parameters analysed two weeks later, i. e. in 8-week-old pigs. In conclusion, chronic exposure to high cortisol levels in HPAhi pigs does not negatively impact immunity.

Exploring the involvement of Tac2 in the mouse hippocampal stress response through gene networking.

Tachykinin 2 (Tac2) is expressed in a number of areas throughout the brain, including the hippocampus. However, knowledge about its function has been only well explored in the hypothalamus in the context of reproductive health. In this study, we identified and validated increased hippocampal Tac2 mRNA expression in response to chronic mild stress in mice. Expression quantitative trait locus (eQTL) analysis showed Tac2 is cis-regulated in the hippocampus. Using a systems genetics approach, we constructed a Tac2 co-expression network to better understand the relationship between Tac2 and the hippocampal stress response. Our network identified 69 total genes associated with Tac2, several of which encode major neuropeptides involved in hippocampal stress signaling as well as critical genes for producing neural plasticity, indicating that Tac2 is involved in these processes. Pathway analysis for the member of Tac2 gene network revealed a strong connection between Tac2 and neuroactive ligand-receptor interaction, calcium signaling pathway, as well as cardiac muscle contraction. In addition, we also identified 46 stress-related phenotypes, specifically fear conditioning response, that were significantly correlated with Tac2 expression. Our results provide evidence for Tac2 as a strong candidate gene who likely plays a role in hippocampal stress processing and neural plasticity.

On association of the lethal yellow (AY) mutation in the agouti gene with the alterations in mouse brain and behavior.

Lethal yellow (AY) mutation causes obesity and type-2 diabetes in mice. Here we studied the effect of the AY mutation on the brain and behavior. The experiments were carried out on adult (11-12 weeks old) males of AY/a mice and their wild-type littermates (a/a). Mice of AY/a and a/a genotypes did not differ in their home cage activity, sleep, food and water consumption, learning ability in the Morris water maze, anxiety in the open field and elevated plus-maze, as well as in the level of monoamines, metabolites and some genes expression in the brain. At the same time, the fat mass, depressive-like immobility in the forced swim and tail suspension tests were significantly increased in AY/a mice compared with a/a ones. Magnetic resonance imaging revealed a significant reduction of cortex volume in AY/a mice. The level of mRNA of Ptpn5 gene encoding striatal enriched tyrosine phosphatase in the frontal cortex of AY/a mice was significantly elevated compared with their wild-type littermates. This is the first report on the alterations in the brain and behavior in the AY/a mouse line. It is tempting to speculate that this mouse line can serve as a new and useful preclinical model to study neurobehavioral complications associated with obesity and type-2 diabetes.

Genetic Contribution to Initial and Progressive Alcohol Intake Among Recombinant Inbred Strains of Mice.

We profiled individual differences in alcohol consumption upon initial exposure and during 5 weeks of voluntary alcohol intake in female mice from 39 BXD recombinant inbred strains and parents using the drinking in the dark (DID) method. In this paradigm, a single bottle of 20% (v/v) alcohol was presented as the sole liquid source for 2 or 4 h starting 3 h into the dark cycle. For 3 consecutive days mice had access to alcohol for 2 h followed by a 4th day of 4 h access and 3 intervening days where alcohol was not offered. We followed this regime for 5 weeks. For most strains, 2 or 4 h alcohol intake increased over the 5-week period, with some strains demonstrating greatly increased intake. There was considerable and heritable genetic variation in alcohol consumption upon initial early and sustained weekly exposure. Two different mapping algorithms were used to identify QTLs associated with alcohol intake and only QTLs detected by both methods were considered further. Multiple suggestive QTLs for alcohol intake on chromosomes (Chrs) 2, 6, and 12 were identified for the first 4 h exposure. Suggestive QTLs for sustained intake during later weeks were identified on Chrs 4 and 8. Thirty high priority candidate genes, including Entpd2, Per3, and Fto were nominated for early and sustained alcohol intake QTLs. In addition, a suggestive QTL on Chr 15 was detected for change in 2 h alcohol intake over the duration of the study and Adcy8 was identified as a strong candidate gene. Bioinformatic analyses revealed that early and sustained alcohol intake is likely driven by genes and pathways involved in signaling, and/or immune and metabolic function, while a combination of epigenetic factors related to alcohol experience and genetic factors likely drives progressive alcohol intake.

Methodologies for Assessing Disease Tolerance in Pigs.

Features of intensive farming can seriously threaten pig homeostasis, well-being and productivity. Disease tolerance of an organism is the adaptive ability in preserving homeostasis and at the same time limiting the detrimental impact that infection can inflict on its health and performance without affecting pathogen burden per se. While disease resistance (DRs) can be assessed measuring appropriately the pathogen burden within the host, the tolerance cannot be quantified easily. Indeed, it requires the assessment of the changes in performance as well as the changes in pathogen burden. In this paper, special attention is given to criteria required to standardize methodologies for assessing disease tolerance (DT) in respect of infectious diseases in pigs. The concept is applied to different areas of expertise and specific examples are given. The basic physiological mechanisms of DT are reviewed. Disease tolerance pathways, genetics of the tolerance-related traits, stress and disease tolerance, and role of metabolic stress in DT are described. In addition, methodologies based on monitoring of growth and reproductive performance, welfare, emotional affective states, sickness behavior for assessment of disease tolerance, and methodologies based on the relationship between environmental challenges and disease tolerance are considered. Automated Precision Livestock Farming technologies available for monitoring performance, health and welfare-related measures in pig farms, and their limitations regarding DT in pigs are also presented. Since defining standardized methodologies for assessing DT is a serious challenge for biologists, animal scientists and veterinarians, this work should contribute to improvement of health, welfare and production in pigs.

Time course study of the response to LPS targeting the pig immune gene networks.

BACKGROUND: Stress is a generic term used to describe non-specific responses of the body to all kinds of challenges. A very large variability in the response can be observed across individuals, depending on numerous conditioning factors like genetics, early influences and life history. As a result, there is a wide range of individual vulnerability and resilience to stress, also called robustness. The importance of robustness-related traits in breeding strategies is increasing progressively towards the production of animals with a high level of production under a wide range of climatic conditions and management systems, together with a lower environmental impact and a high level of animal welfare. The present study aims at describing blood transcriptomic, hormonal, and metabolic responses of pigs to a systemic challenge using lipopolysaccharide (LPS). The objective is to analyze the individual variation of the biological responses in relation to the activity of the HPA axis measured by the levels of plasma cortisol after LPS and ACTH in 120 juvenile Large White (LW) pigs. The kinetics of the response was measured with biological variables and whole blood gene expression at 4 time points. A multilevel statistical analysis was used to take into account the longitudinal aspect of the data. RESULTS: Cortisol level reaches its peak 4 h after LPS injection. The characteristic changes of white blood cell count to LPS were observed, with a decrease of total count, maximal at t=+4 h, and the mirror changes in the respective proportions of lymphocytes and granulocytes. The lymphocytes / granulocytes ratio was maximal at t=+1 h. An integrative statistical approach was used and provided a set of candidate genes for kinetic studies and ongoing complementary studies focused on the LPS-stimulated inflammatory response. CONCLUSIONS: The present study demonstrates the specific biomarkers indicative of an inflammation in swine. Furthermore, these stress responses persist for prolonged periods of time and at significant expression levels, making them good candidate markers for evaluating the efficacy of anti-inflammatory drugs.

Differentially expressed genes and gene networks involved in pig ovarian follicular atresia.

Ovarian folliculogenesis corresponds to the development of follicles leading to either ovulation or degeneration, this latter process being called atresia. Even if atresia involves apoptosis, its mechanism is not well understood. The objective of this study was to analyze global gene expression in pig granulosa cells of ovarian follicles during atresia. The transcriptome analysis was performed on a 9,216 cDNA microarray to identify gene networks and candidate genes involved in pig ovarian follicular atresia. We found 1,684 significantly regulated genes to be differentially regulated between small healthy follicles and small atretic follicles. Among them, 287 genes had a fold-change higher than two between the two follicle groups. Eleven genes (DKK3, GADD45A, CAMTA2, CCDC80, DAPK2, ECSIT, MSMB, NUPR1, RUNX2, SAMD4A, and ZNF628) having a fold-change higher than five between groups could likely serve as markers of follicular atresia. Moreover, automatic confrontation of deregulated genes with literature data highlighted 93 genes as regulatory candidates of pig granulosa cell atresia. Among these genes known to be inhibitors of apoptosis, stimulators of apoptosis, or tumor suppressors INHBB, HNF4, CLU, different interleukins (IL5, IL24), TNF-associated receptor (TNFR1), and cytochrome-c oxidase (COX) were suggested as playing an important role in porcine atresia. The present study also enlists key upstream regulators in follicle atresia based on our results and on a literature review. The novel gene candidates and gene networks identified in the current study lead to a better understanding of the molecular regulation of ovarian follicular atresia.

Acute social stress-induced immunomodulation in pigs high and low responders to ACTH.

Pig husbandry is known as an intensive breeding system, piglets being submitted to multiple stressful events such as early weaning, successive mixing, crowding and shipping. These stressors are thought to impair immune defences and might contribute, at least partly, to the prophylactic use of antibiotics. Robustness was recently defined as the ability of an individual to express a high-production potential in a wide variety of environmental conditions. Increasing robustness thus appears as a valuable option to improve resilience to stressors and could be obtained by selecting piglets upon their adrenocortical activity. In this study, we aimed at depicting the consequences of an acute social stress on the immune capacity of piglets genetically selected upon divergent hypothalamic-pituitary-adrenocortical (HPA) axis activity. For this purpose, we monitored neuroendocrine and immune parameters, in high- (HPAhi) and low- (HPAlo) responders to ACTH, just before and immediately after a one-hour mixing with unfamiliar conspecifics. As expected, stressed piglets displayed higher levels of circulating cortisol and norepinephrine. Blood cell count analysis combined to flow cytometry revealed a stress-induced leukocyte mobilization in the bloodstream with a specific recruitment of CD8α+ lymphocytes. Besides, one-hour mixing decreased LPS-induced IL-8 and TNFα secretions in whole-blood assays (WBA) and reduced mononuclear cell phagocytosis. Altogether, our data demonstrate that acute social stress alters immune competence of piglets from both groups, and bring new insights in favour of good farming practices. While for most parameters high- and low-responders to ACTH behaved similarly, HPAhi piglets displayed higher number of CD4+ CD8α- T cells, as well as increased cytokine production in WBA (LPS-induced TNFα and PIL-induced IL-8), which could confer them increased resistance to pathogens. Finally, a principal component analysis including all parameters highlighted that overall stress effects were less pronounced on piglets with a strong HPA axis. Thus, selection upon adrenocortical axis activity seems to reduce the magnitude of response to stress and appears as a good tool to increase piglet robustness.

Time course of the response to ACTH in pig: biological and transcriptomic study.

BACKGROUND: HPA axis plays a major role in physiological homeostasis. It is also involved in stress and adaptive response to the environment. In farm animals in general and specifically in pigs, breeding strategies have highly favored production traits such as lean growth rate, feed efficiency and prolificacy at the cost of robustness. On the hypothesis that the HPA axis could contribute to the trade-off between robustness and production traits, we have designed this experiment to explore individual variation in the biological response to the main stress hormone, cortisol, in pigs. We used ACTH injections to trigger production of cortisol in 120 juvenile Large White (LW) pigs from 28 litters and the kinetics of the response was measured with biological variables and whole blood gene expression at 4 time points. A multilevel statistical analysis was used to take into account the longitudinal aspect of the data. RESULTS: Cortisol level reached its peak 1 h after ACTH injection. White blood cell composition was modified with a decrease of lymphocytes and monocytes and an increase of granulocytes (F D R<0.05). Basal level of cortisol was correlated with birth and weaning weights. Microarray analysis identified 65 unique genes of which expression responded to the injection of ACTH (adjusted P<0.05). These genes were classified into 4 clusters with distinctive kinetics in response to ACTH injection. The first cluster identified genes strongly correlated to cortisol and previously reported as being regulated by glucocorticoids. In particular, DDIT4, DUSP1, FKBP5, IL7R, NFKBIA, PER1, RGS2 and RHOB were shown to be connected to each other by the glucocorticoid receptor NR3C1. Most of the differentially expressed genes that encode transcription factors have not been described yet as being important in transcription networks involved in stress response. Their co-expression may mean co-regulation and they could thus provide new patterns of biomarkers of the individual sensitivity to cortisol. CONCLUSIONS: We identified 65 genes as biological markers of HPA axis activation at the gene expression level. These genes might be candidates for a better understanding of the molecular mechanisms of the stress response.

Transcriptome profiling of sheep granulosa cells and oocytes during early follicular development obtained by laser capture microdissection.

BACKGROUND: Successful achievement of early folliculogenesis is crucial for female reproductive function. The process is finely regulated by cell-cell interactions and by the coordinated expression of genes in both the oocyte and in granulosa cells. Despite many studies, little is known about the cell-specific gene expression driving early folliculogenesis. The very small size of these follicles and the mixture of types of follicles within the developing ovary make the experimental study of isolated follicular components very difficult.The recently developed laser capture microdissection (LCM) technique coupled with microarray experiments is a promising way to address the molecular profile of pure cell populations. However, one main challenge was to preserve the RNA quality during the isolation of single cells or groups of cells and also to obtain sufficient amounts of RNA.Using a new LCM method, we describe here the separate expression profiles of oocytes and follicular cells during the first stages of sheep folliculogenesis. RESULTS: We developed a new tissue fixation protocol ensuring efficient single cell capture and RNA integrity during the microdissection procedure. Enrichment in specific cell types was controlled by qRT-PCR analysis of known genes: six oocyte-specific genes (SOHLH2, MAEL, MATER, VASA, GDF9, BMP15) and three granulosa cell-specific genes (KL, GATA4, AMH).A global gene expression profile for each follicular compartment during early developmental stages was identified here for the first time, using a bovine Affymetrix chip. Most notably, the granulosa cell dataset is unique to date. The comparison of oocyte vs. follicular cell transcriptomes revealed 1050 transcripts specific to the granulosa cell and 759 specific to the oocyte.Functional analyses allowed the characterization of the three main cellular events involved in early folliculogenesis and confirmed the relevance and potential of LCM-derived RNA. CONCLUSIONS: The ovary is a complex mixture of different cell types. Distinct cell populations need therefore to be analyzed for a better understanding of their potential interactions. LCM and microarray analysis allowed us to identify novel gene expression patterns in follicular cells at different stages and in oocyte populations.

Molecular genetics of hypothalamic-pituitary-adrenal axis activity and function.

The hypothalamic-pituitary-adrenocortical (HPA) axis is a major neuroendocrine system involved in the regulation of numerous physiological processes and in adaptation to stress. A wide range of variability can be observed in all the components of the system, and the contribution of genetic factors has been shown in the central regulation of the axis, the production of glucocorticoid hormones by the adrenal cortex, their bioavailability, and the efficiency of their action at the level of receptor and postreceptor mechanisms. Numerous molecular polymorphisms have been described that contribute to physiological variation as well as to HPA axis-related pathological conditions. Although most studies focus on single gene polymorphisms, future studies should aim to integrate the different sources of variation into a systems genetic model to take into account the strong interdependence of the different components of the axis.